DIFFUSION IN RANDOM ENVIRONMENT AND 
THE RENEWAL THEOREM 



DIMITRIOS CHELIOTIS 



Abstract. According to a theorem of S. Schumacher and T. Brox, for a dif- 
fusion X in a Brownian environment it holds that (Xt — b\ ogt )/log 2 t — > 
in probability, as t — > oo, where b is a stochastic process having an explicit 
description and depending only on the environment. We compute the distri- 
bution of the number of sign changes for b on an interval [1, x] and study some 
of the consequences of the computation; in particular we get the probability of 
b keeping the same sign on that interval. These results have been announced 
in 1999 in a non-rigorous paper by P. Le Doussal, C. Monthus, and D. Fisher 
and were treated with a Renormalization Group analysis. We prove that this 
analysis can be made rigorous using a path decomposition for the Brownian 
environment and renewal theory. Finally, we comment on the information 
these results give about the behavior of the diffusion. 



On the space W := C(R) consider the topology of uniform convergence on com- 
pact sets, the corresponding er-field of the Borel sets, and P the measure on W under 
which the coordinate the processes {w(t) : t > 0}, {w(—t) : t > 0} are independent 
standard Brownian motions. 

Also let ft := C([0, +oo)), and equip it with the cr-ficld of Borel sets derived 
from the topology of uniform convergence on compact sets. For w € W, we denote 
by P w the probability measure on f2 such that : t > 0} is a diffusion with 

w(0) = and generator 



The construction of such a diffusion is done with scale and time transformation from 
a one-dimensional Brownian motion (see, e.g., |14j . |16j). Using this construction, 
it is easy to see that for P-almost all w s W the diffusion does not explode in finite 
time; and on the same set of u>'s it satisfies the formal SDE 



where (3 is a one-dimensional standard Brownian motion. 
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1. Introduction 




(1) 



dw{t) = df3(t) - \w'{uj{t))dt, 
w(0) = 0, 
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Then consider the space W x fi, equip it with the product er-field, and take the 
probability measure defined by 

dP(w,u) = dP w (uj)dP{w). 

The marginal of V in gives a process that is known as diffusion in a random 
environment; the environment being the function w. 
S. Schumacher f |14M13| 1 proved the following result. 

Fact 1. There is a process b : [0, oo) x W — > R such that for the formal solution tu 
of Q it holds 

(2) - - N9 - hiw^V) -> in V as t -> +oo, 

(logty 

where for r > we let «/ r )(s) = r~ 1 w(sr 2 ) /or a/Z s £ M. 

We will define the process & soon. This result shows the dominant effect of the 
environment, through the process b, on the asymptotic behavior of the diffusion. 
The results we prove in this paper concern the process b. In subsection 11.21 we 
commend on their implications for the behavior of the diffusion itself. 

Besides this diffusion model, there is a discrete time and space analog, known 
as Sinai's walk, which was studied first. Sinai's pioneering paper |17| identified the 
role of the process b in the analogous to © limit theorem for the walk. Then S. 
Schumacher proved in ^31 ( see also [13] f° r * ne results without the proofs) a more 
general statement than the above proposition where the environment w was not 
necessarily a two sided Brownian motion, while T. Brox pQ gave a different proof 
in the Brownian case. H. Kesten [Hj computed the density of b\ in the case we 
consider, and Tanaka |19j generalized the computation to the case that w is a two 
sided symmetric stable process. Localization results have been given for the Sinai 
walk by Golosov (0, actually, for the reflected walk) and for the diffusion model 
by Tanaka Also Tanaka (JH], El) studied the cases where the environment is 
non-positive reflecting Brownian motion, non-negative reflecting Brownian motion, 
or Brownian motion with drift. Finer results on the asymptotics of Sinai's have been 
obtained by Z. Shi and Y. Hu. A survey of some of them as well as a connection 
between Sinai's walk and diffusion in random environment is given in |16| . Another 
connection is established in |15j . 

In ^U], P. Le Dousal, C. Monthus, and D. Fisher proposed a new method for 
tackling questions related to asymptotic properties of Sinai's walk, and using it 
they gave a host of results. The method is a Rcnormalization Group analysis and 
it has consequences agreeing with rigorously proved results (e.g., 0, |H]) ■ This 
is the starting point of the present paper. In the context of diffusion in random 
environment, we show how one can justify the Rcnormalization Group method using 
two tools. The first is a path decomposition for a two sided standard Brownian 
motion; the second is the renewal theorem. Our main results illustrate the use of 
the method and the way we justify it. 

The structure of the paper is as follows. In the remaining of the introduction we 
state our results. In Section[2|wc provide all the necessary machinery for the proofs, 
which arc given in Section [3J Some technical lemmata that we use are proved in 
Section HI 

We begin by defining the process b. 
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For a function w : R — ► R, x > and i/o £ R we say that w admits an x- 
minimum at yo if there are a,|3 € I with a < yo < (3, w(yo) = m{{w(y) : y G 
[a, 0]} and w(a) > w(yo) + x, w(j3) > w(yo) + x. We say that w admits an 
x- maximum at yo if — w admits an x- minimum at yo. 



a yo 13 
Figure 1. w admits an x-minimum at yo- 

For convenience, we will call a point where w admits an x-maximum or x- 
minimum an x-maximum or an x-minimum respectively. 
We denote by R x (w) the set of x-extrema of w and define 

{For every x > the set R x (w) has no accumulation point in R, 
w G W : it is unbounded above and below, 

and the points of x-maxima and x-minima alternate. 

Thus, for w G Wi and x > we can write R x (w) — {xh{w,x) : k G Z} with 
(x k (w,x)) keZ strictly increasing, x (w,x) < < x\(w,x) , lim^-oo x k (w, x) = 
— oo, limfc^oo Xk{w,x) = oo. It holds that P(Wi) = 1, and the easy proof of this 
fact is given in Lemma El 

Definition 1. The process b : [0, +oo) x W — > R is defined for x > and w G Wi 

as 



Xo(w, x) if x a (w, x) is an x- minimum, 
xi(w, x) else, 



and b x {w) = if x = or to G W \ Wi- 

Remark 1. In the definition of b x (w) wc do not make use of the entire sequence 
of x-extrema. The reason we introduce this sequence is that we plan to study the 
evolution of the process b as x increases. Since R x (w) C R x (w) for x < x, the later 
values of b.(w) are elements of R x (w). For x large enough, the points xq(w,x), 
xi(w,x) will not be x-extrema. 

Remark 2. We will decompose the process uu at the endpoints of the intervals 
{[xfc(w, x), Xfc+i (w, x)] : k G Z} and study its restriction to each of them. Of course 
[xo(w, x), Xi (w, x)] has a particular importance for the process 6; and it is in the 
study of w\ [xo(w, x), x\ (w, x)] that the renewal theorem enters (see Lemmanj. 

Remark 3. It is clear that b satisfies b x (—w(c-)) — -b ax (w) for all a, c, x > 0, and 
w G W. So that the quantity &i(ui^° st ^) appearing in @ equals also bi og t(w) / (\ogt) 2 
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1.1. Sign changes of b. For x > 1 define on Wi the random variable 

fc(x) = f/= times b.(w) has changed sign in [l,x]. 
The main result of the paper is the computation of the generating function of k(x). 
Theorem. For x > 1 and z e C wrei/i |z| < 1, it Zio/ds 

(3) E(z fc(a; )) = ci(z) :c Al « + c 2 (z) a^W, 
w/iere 

. , , -3 + V5 + 4z x . , -3 - V5 + 4z 
A i(2) = g ' Aa ( z ) = 2 ' 

and 

Cl (z) = ((z-l)/3-A 2 (z))/(A 1 (z)-A 2 (z)), c 2 (z) = (-(z-l)/3+A 1 (z))/(A 1 (z)-A 2 (z)). 

From this we extract several corollaries. Corollary Q] and Corollary [3] are imme- 
diate, while the rest require some work and are proved in section [5J 

Corollary follows by taking z — > in 

Corollary 1. 

(4) F(b.(w) keeps the same sign in [1, x])/x ( " 3+v/5)/2 -> 1/2 + 7V5/30 
as x — > +oo. 

Corollary 2. TTie increasing process of points {Xk)k>\ where b changes sign in 
[1, +oo) has the form X^ = X±n ■ ■ -rfc-i, k > 2, where X\ is the smallest such 
point and the rj 's are i.i.d wrf/i density 

(5) /( r ) = 7 ^_( r Ai-l_ r A a -l) ) r > 1, 

Ai — A 2 

where Ai = Ai(0), A 2 = A 2 (0). 
Now observe that 

k(t) = sup{?i G N : X n < t} = sup{n G N : log Xi + logn H h logr„_i < logi}. 

Since E(logri) = 3 and logXi is finite a.s. (e.g., by Corollary^ E(logXi) < +oo), 
the next statement follows from renewal theory. 

Corollary 3. fc(t)/logf -> 1/3 as t ^ +co F-a.s. 

Corollary |21 allows the following strengthening of the above theorem. 

Corollary 4. Relation © holds for all z G C \ (-co, -5/4]. 

For t G (0, oo) consider the random variable k(e t )/t and let fit be its distribution 
measure. Then the following holds. 

Corollary 5. The family of measures (/it)t>0 satisfies a Large Deviation Principle 
with speed t and good rate function 



I{x) 



x log (2x(x + y/x 2 + 5/4)) + § - (x + y/x 2 + 5/4) if x G [0, +oo), 
+oo if x £ (—oo, 0). 
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In ^Hj, Corollary |21 appears in paragraph IV. B with a different justification. We 
state it here because we need it for the proof of Corollary The large deviation 
result of Corollary [S] is the precise mathematical interpretation of the discussion in 
paragraph IV. A of the same paper. 

In Corollary ^ the condition k(x) = means that the process b tends to keep 
the diffusion away form zero in the time interval [e, e x ] (since the diffusion localizes 
around b, and b keeps sign on [1, x]). For the event that the diffusion hits zero there 
are two interesting relevant papers. The first one is by Y. Hu who treats the 
annealed asymptotics of the first time of hitting zero after time t as t — > +oo . The 
second is by F. Comets and S. Popov and refers to a related model. That is, it 
considers a process {X t )t>o on Z that runs in continuous time in an environment 
us satisfying the conditions of the Sinai model and studies, among other things, the 
asymptotics of the quenched probability P w (X t = 0\Xq = 0) as t — ► +oo. 

1.2. The process b and the diffusion. The results of the previous subsection 
concern the process &, which is a functional of the environment. And our motivation 
for studying b was the localization results involving this process (the simplest being 
Fact Q an d keep in mind Remark yj. An obvious question is what we can infer 
about the behavior of the diffusion from our results. 

Using the representation of the diffusion as a time and scale change of Brownian 
motion, one can show easily that the diffusion is recurrent and is a regular point 
for (0, +oo) and (— oo, 0). Consequently, for all c > 0, the diffusion visits in the 
time interval [c, +oo) and in its first visit there it scores an infinite number of sign 
changes. So there can be no direct connection with the corresponding number for 
the process b. One can consider, say, the number of sign changes of the diffusion 
between times where the diffusion achieves a positive or negative record value. This 
number is finite on compact intervals of (0, +oo) but still not related with the sign 
changes of b (and it is not hard to sec this). 

When &i og . jumps to a new value, what happens is not that just the diffusion goes 
through that value shortly before or after that. As it is known (see, e.g., |2()j). 
the impact of the jump is that the diffusion goes to that value and it is trapped in 
its neighborhood for a large amount of time. The way to detect the approximate 
location of such values when we observe the diffusion in real time (i.e., at time t we 
know w|[0, i]) is to find the site the diffusion has spent the most time thus far. 

To make the last statement precise, we use the local time process {L u (t, x) :t> 
0, x € R} that corresponds to the diffusion ui. This process is jointly continuous, 
and with probability one satisfies 



for all £ > and any bounded Borcl function / S R . 

For a fixed t > 0, the set $(t) := {x G R : L u {t,x) = sup eK L u (t, y)} of the 
points with the most local time at time t is nonempty and compact. Any point 
there is called a favorite point of the diffusion at time t. Define F : (0, +oo) — > R 
with F(t) = inf^(i), the smallest favorite point at time t. Also, for x > 1 define 
the interval I(x) := (x — (\ogx) 5 ,x + (logx) 5 ). 

In a work in progress we expect to prove the following. 
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Claim. With V probability one, there is a strictly increasing sequence (t n (u), w)) n >i 
converging to infinity so that if we denote by (x n {w,ti)) n >i the sequence of con- 
secutive values that b\ og . takes on the interval (ti,+oo) (remember that b is a step 
function in any interval [x, +00) with x>0), then 



and x n = fciogt for some t G (t n ,t n +i)- We abbreviated t n (u),w),x n (w,ti) tot n ,x n . 

Observe that for big x, the interval I(x) is a relatively small neighborhood of 
x. Thus, the claim says that, after some point, the function F "almost tracks" 
the values of the process b\ og . with the same order and at about the same time. 
Consequently, the number of sign changes of 6i og . on an interval (s,t) (with s,t 
large) would correspond to the number of sign changes of F on approximately the 
same interval. Or, more precisely, the number of sign changes of F on (s, t) and 
the corresponding number for b\ og . differ by at most two. It is easy to see that the 
last statement follows from the Claim above. 



As a first step towards the study of the process b, we look at the law of the 
Brownian path between two consecutive x-extrema as well as the way these pieces 
are put together to constitute the entire path. 

The first piece of information is provided by Proposition of §1 in |12j . 

Fact 2. For every x > 0, the times of x-extrema of a Brownian motion (wt ■ t G 
K), wo — build a stationary renewal process R x {w) = {xk : k G Z} with (xk)kez 
strictly increasing and Xq < < X\. The trajectories between consecutive x-extrema 
(w Xk +t — Wx k '■ t S [0, Xk+i — Xk]), k € Z are independent and for k ^ identically 
distributed (up to changes of sign). 

In Lemma of §1 of the same paper |12| a description of each such trajectory is 
given which we quote (see Figure |2J). 



F((t n ,t n+1 )) C I(x n ) for all n > 1, 



2. Preliminaries 



w 




s 



Figure 2. The graph of w until M — w hits x. 



For x, t > let 



a 



T 



sup{w s : s G [0, <]}, 
min{f : M t — Wt + x} 

max{s G [0, r] : w s = 
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Fact 3. The following hold: 

(i) The two trajectories (wt '■ t € [0, a]) and ((3 — w a +t '■ t € [0,t — a]) are 
independent. 

(ii) (3 has exponential distribution with mean x. 

(iii) The Laplace transform of the law of a given (3 is 

E[e- sa \/3 = y] = exp{-(y / 'x)<j){sx 2 )} , s > 0. 

(iv) The Laplace transform of the law of t — o is 

E [ e -s(r-<x)] = ^( sa ,2) , S > 0, 

where 4>(s) = \/2s coth ^/2s — 1 and i()(s) = V2s/ sinh\/2S' 

We call the translation (w — w(xk))\[xk, Xk+i] of the trajectory of w between two 
consecutive rr-extrema an 2>slope (or a slope, when the value of x is clear or irrel- 
evant), a slope that takes only non-negative values an upward slope, and a slope 
taking only non-positive values a downward slope. We call (w — w(xq))\ [xq, x\] 
the central x-slope. 

In the following we will use the operation of "gluing together" functions defined 
on compact intervals. For two functions / : [a, 0] — > M, g : [7, 5] — > K, by gluing g 
to the right of / we mean that we define a new function j : [a, f3 + 5 — 7] — > R with 

■( ) _ ff( x ) for x <E [a,/3], 

' \f(J3)+ g(x-P + 7) -9(7) facxe\p,i3 + 6-i\. 

It is clear that for all k 7^ 0, if (w — w(xk))\[xk, Xk+i] is an upward x-slope, then 
it is obtained by gluing together a trajectory of type (w t : t G [0, a]) to the right of 
a trajectory of type (/3 — w a+t : t € [0, r — a]) and then translating the resulting 
path so that it starts at (x^, 0) . Similarly for a downward slope. 

For any x-slope T : [a, 0\ -> E we call i(T) := /3 - a , fo(T) := \T(j3) - T(a)\ the 
length and the height of the slope respectively, and we denote by r](T) := h(T)—x 
the "excess height" of T. Also we denote by |T|, 9(T), the slopes with domains 
[a, p], [0,13 -a] and values |T|(-) = |T(-)|, 6>(T)(-) = T(a + •) respectively. 

For any slope T, the slope \6(T)\ is in the set S defined by 

{/ is a function such that there is a 1(f) > with "1 
/ C [0, +00) 2 : Domain(/) = [0, 1(f)], f continuous, I . 

and = /(0) < f(x) < /(!(/)) Vx € [0, /(/)] J 

On 5 we define a topology for which the base of neighborhoods of an element f E S 
is the collection of all sets of the form 

{geS: \l(g) - l(f)\ < e and \f(t 1(f)) - g(t l(g))\ < e for all t e [0, 1]}. 

With this topology, S is a Polish space. Equip S with the Borel cr-algcbra, and de- 
fine the measures m r x , m% the first to be the distribution of 9(\(w — w(x 1 ))\\[xi,X2]) 
and the second to be the distribution of 6(\(w — w(xq)) \ | [xq, x\]) (the superscripts 
r and c standing for renewal and central). 

In the remaining part of this section we compute the distribution of the length 
and excess height of a slope picked from m r x or m c x . We assume x = 1 since the 
scaling property of Brownian motion gives the corresponding results for the case 
x ^ 1. 
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Let T be a slope picked from m\. Earlier we described the way an upward slope 
is formed. Its length is the sum of two independent independent random variables 
Zi, Zi with Z\ lc = a, Z2 l< = t — a (where we take x = 1 in their definitions just 
before Fact OJl ■ But from (i) of Fact er and r — a are independent. Thus, 

l(T) t. By definition, r is the time where the reflected process M — w hits one. 
This reflected process has the same law as \w\, and the Laplace transform of the 
time it first hits one is known as 

(6) E[e- Xl ^} = (cosh^)- 1 for A > 0. 

Also K(l(T)) = E(t) = 1. Using the Laplace inversion formula (see PS 531) 
we find that the density of I = l(T) is 

(7) / i (x) = |^(-l) fc (fc + i)cx P h^(fc + i) 2 x] ,x>0. 

feez 

We note for future reference that, by (ii) of Fact [3 for any a > 0, the excess height 
of a slope picked from m r a is exponential with mean a; i.e., it has density 

(8) Pa (x) = a - 1 e- x/a , x > 0. 

Now let To be a slope picked from mf. More specifically, take To to be the cen- 
tral 1-slope. Observe that, by virtue of Fact [3 we can "start a renewal process 
at —00" with i.i.d. alternating upward and downward slopes, and ask what the 
characteristics are of the slope covering zero. The renewal theorem says that the 
length of the slope covering zero is picked from the distribution of I given in l(7)l 
with size-biased sampling. Once the length, say z, is picked, we expect that the 
remaining characteristics of the slope, ignoring direction (upward or downward), 
are determined by the law of T|Z(T) = z under m\. We give a formal proof of this. 
Notice that a regular conditional distribution for the random variable T given the 
(7-ficld a(l(T)) exists because the space iS is Polish. 

Lemma 1. For any measurable subset A of S, it holds 

00 

(9) n\d(T )\ e A) = J P(T e A I l(T) = z)zf l (z) dz, 



where T has under P distribution m\. 

Proof. Let Fi be the distribution function of / = l(T), and for t € K let T(t) be the 
1-slopc around t. That is, the slope whose domain contains t. Then P(|0(To)| G 
A) = ¥(\6(T(t))\ e A) for all t > because 9(T ) is the same as the image under 6 of 
the slope around t for [w s -t— W—t : s € R), and the latter process is again a standard 
two sided Brownian motion. Now let (Y n ) n >o be an independent sequence of slopes 

with Y n l =" ( — 1)™ +1 T. Glue them sequentially to get a function / in C([0, +00)) 
with /(0) = 0, and denote by T(t) the slope around t, for t > 0. 

If we take a as defined just before Fact [3] with x = 1 in the definition of r there, 
then it holds 

P(|0(T o )| e A) = n\0(T(t))\ e A) = F(\0(T(t))\ e A, a < t)+F(\9(T(t))\ €A,a>t) 
and 

n\9(T(t))\ eA,cj<t)= f ¥(\9(f(t -y))\e A)f a (y) dy, 

Jo 
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where f a is the density of a. We will take t —* +00 and finish with the proof by 
showing that the limit lim t _> +oc P(\9(T(t))\ <G A) exists and 

00 

lim P(|0(T (t))\ G A) = [p(T G A I l(T) = z)zfAz) dz. 

t^+OQ J 



To see this, define g(t) := P(|0(T(i))| G A) for t > 0. Then 

g(t)=F(T€A, l(T)>t) + j ¥{\6{f (t - s))\ € A) dFi(s) 



/ 



»(Te A l(T)>t)+ / fl(t-«) dF;(s). 





The distribution of Z is nonarithmctic with mean value 1. By the renewal theorem 
([Sj, Chapter 3, statement (4.9)), it follows that the linii_> +00 g(t) exists and 

t lim g(t) = J P(T e A, l(T) > s) ds = J J P(T G A | ?(T) = z)/;(z) dzds 

s 

00 z 00 

?(T eA \ l(T) = z)fi(z) dsdz = Jp{T eA \ l(T) = z)zfi(z) dz. 
00 



Now we apply Lemma ^ to obtain the distribution of the length and height of 
the central 1-slope. 

• Regarding the length of To, observe that for x > the set A := {T G S : 
l(T) < x} is open and 

00 X 

P(/(T ) < x) = J P(l(T) < x I l(T) = z)zf l (z) dz = J z f,(z) dz. 

So that /(To) has the density 

(10) f l{To) (x) = xft(x) toxx>0, 

which is the size-biased sampling formula from renewal theory. 

• Regarding the excess height of To, observe that for x > the set A := {T G 
S : rj(T) < x} is open and 

00 00 

P(r/(T ) <x) = J P(r?(T) < x \ l(T) = z)zfi(z) dz = J P(r?(T) < x, l(T) = z)z dz. 


Differentiating with respect to x we get the density of ?7(Tb) as 

B (?7(T) = x, l(T) = z)z dz = E (l(T) | r?(T) = x) e' x . 





10 



DIMITRIOS CHELIOTIS 



From (iii), (iv) of Fact El 

E(e~" (T) I V(T) = x) = i>(t) e- x ^\ 

After some calculations, JjV'W e~ x ^^\ t= o = — (2x + l)/3. The derivative can 
move inside the expectation on the left hand side of the above relation, giving 
— E(Z(T) | r](T) = x), due to the monotone convergence theorem; which applies 
because l(T) > and the function (x i— > (1 — e~ ax )/x) is nonnegative and decreasing 
in R for any a > 0. Therefore the density of 77 (To) is 

(11) W*) ^ 2 *^ 6 "* fo ra; >o. 

The last bit of information needed to achieve our goal, stated at the first sentence 
of this section, is the direction of the central 1-slopc (i.e., upward or downward) 
and its location with respect to zero. By symmetry, To is an upward slope with 
probability 1/2, and from exercise 3.4.7 of it follows easily that given the length 
I of the slope (w ~ w(xq))\[xo, x±] around zero, the distance of zero from xq is 
uniformly distributed in [0,1]. 

3. Proof of the Theorem and the Corollaries 

For x > and w £ Wi with set of x-extrema R x (w) = {xk ■ k £ Z} define 

A x (w) := {(w - w(xk)j\[xk,Xk+i] ■ k £ zj. 

We refer to the parameter x as time since we are going to study the evolution of 
A x {w) as x increases. A x {w) is the set of slopes at time x. Roughly, as x increases, 
the slopes that have height smaller than x are absorbed into greater ones. 
For any Lebcsgue measurable set S C [0, +00) , x > 1, and k £ N we define 

In A x {w) the central slope has excess height y £ S 
and b. (w) has changed sign k times in [l,x] 

It is important to note that the values of b up to time x are "encoded" in the central 
slope of A x (w). So the number of sign changes of b. (w) in [1, x] can be inferred by 
that slope. 

For x, k fixed, U is a measure that satisfies 

U(x, S, k) < P(ln A x {w) the central slope has excess height y £ S), 

and the right-hand side, considered as a function of S, is a measure absolutely 
continuous with respect to the Lebesgue measure with density {\/x)f n (T ){- /x), 
where f v ,T ) is given in (|ll|l . Therefore the measure on the left-hand side has a 
density also (an element of T 1 ([0, +00))), call it u(x, y, k), and for Lebesgue almost 
all y it holds 

+00 ^ 

( 12 ) "(^ y ' k "> = x^ {To) ^ ~x)' 

Define U : [1, +00) x [0, +00) x N — > [0, +00) with U(x, y, k) := U(x, [y, +00), k). 
U is continuous as is proved in Lemma [7| in section We plan to establish a PDE 
for U. To do this, we look at A x (w) and try to predict how A x+e (w) should look 
like. 

In the transition from A x (w) to A x + E (w), a slope around a point remains the 
same if the excess height of this and the two neighboring slopes are greater than 



U(x,S, k) 
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X 



Figure 3. The decomposition of a piece of the Brownian path in 
x-slopcs. The dots mark the points of cc-extrema. The length x is 
shown on the side. 

e. In case one slope has excess height in [0, e), it does not appear in A x+E {w). For 
example, in Figure|3]the slope To has a height say x + ho with ho € [0,e), and the 
two neighboring slopes T_i, T\, have height greater than x + e. Assume that T_x, 
Ti, have heights x + V\, x + v-i respectively. In A x+£ (w) we know that the slopes 
T_ i , To , T\ will merge to constitute a new slope with height X + V1+X + V2 — x — ho = 
x + v\ +i>2 — ho', i.e., with excess height V\ + V2 — ho — e. The slopes T_2, T2 can stay 
as they arc in the transition from A x (w) to A x+£ or they can be extended if some of 
T_3, T3 has excess height in [0, e). In any case, they don't interfere with T_i, To, T\. 
This simple observation, combined with the renewal structure of the sets A x (w), 
is the basis for the next lemma, which is the first step towards establishing a PDE 
that U solves. We denote by d y U the y derivative of U(x,y,k), and recall that 
p x (v) was defined in JSJ as the density of an exponential with mean x. 

Lemma 2. For x > 1, y > 0, e > 0, k G N we have 



(13) U(x + e, y, k) = U(x, y + e, k)e~ 2e / x + — / p x (v) U(x, (y - v)+,k) dv- 

x Jo 

l{ fc >i}£ d y U(x, Q,k- l)(y/x + l) e -^ x + o(e), 



where for k > 1 we assume that U(x, y, k — 1) is differ entiable in y with continuous 
derivative. The term o(e) depends on x,y,k. 

Proof. The left-hand side of the equation is the probability of an event referring 
to A x+£ (w), and we express it in terms of probabilities referring to A x (w). In 
A x {w) we focus our attention on the seven rr-slopes closest to zero. Denote them 
by Ti, i = — 3, . . . , 3 in the order they appear in the path of w from left to right, 
with T being the central one. The slopes |#(T)|, i = —3, ... ,3 arc independent 
having for i 7^ law m x , and |#(T )| having law m x . The probability of the event 
that at least two of them have excess height in [0,e) is bounded by 21e 2 /x 2 , and 
this is accounted for in the o(e) term in 1)13(1. In the complement of this event, the 
event whose probability appears in the left-hand side of l|13f) happens if and only 
if in A x (w) one of the following three holds (considering what slope, if any, among 
the seven has excess height in [0,e)). 

• At most one of T_3, T_2, T2, T3 has excess height in [0, e), the excess height 
of each of T_i, To, T\ is at least e, and b has changed sign k times in [l,x]. 
In this case the slope around zero is the same for both A x (w) , A x+e (w) . 

• Exactly one of T_i,Tt has excess height in [0,e), T has excess height at 
least e, and b has changed sign k times in [l,x]. In this case b x+e has the 
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same sign as b x . For example, assume that Xi has excess height in [0,e). 
If b x > then b x + £ > b x , while if b x < then b x + £ = b x and simply the 
central slope in A x+£ (w) results from merging Tq,T\,T 2 . 
• k > 1, To has excess height in [0, e), and b has changed sign k — 1 times 
in In this case b x + £ b x < because in A x+£ (w) the central slope 

results from merging T_i,To,Ti and it has different direction than To (i.e., 
if for example To is an upward slope, then the central slope in A x+£ (w) is 
a downward slope). 

The first case has probability 

U(x, y + e 7 k) P(?y_i e [e, +oo)) P(??i e [e, +oo)) = U(x, y + e, k) exp(-2e/a;) 

with7y_i (resp. rji) denoting the excess height of T_i (resp. T\). This follows by the 
independence mentioned above and by the fact that r\-\ and r\\ have exponential 
distribution with mean x, and it expresses the demand that both of the x-slopes 
neighboring the central x-slope have excess height greater than e. In this case, the 
excess heights of T_3, T_2, T2, T3 do not matter. They can't influence the central 
slope in A x+£ (w). 
The second case has probability 

/•S r+00 

2 J J Px{vi)p x (v 2 )U(x,eV (y + s-v 2 +v 1 ) 7 k)dv 2 dv 1 

because, say in the event that T\ has excess height in [0, e), the central slope in 
A x+£ (w) comes from merging Tq,Ti,T 2 of A x (w). Assume that they have excess 
heights u,v\,v 2 respectively. Then the central slope in A x+£ (w) will have excess 
height u — v\ + v 2 ~ e, and the requirement that this is greater than y translates to 
u being greater than (y + v\ — v 2 + e) + . And of course, by assumption, To and T 2 
have excess height greater than e. 

U is continuous as it is proved in Lemma Thus, dividing with e the previous 
integral and taking e — > we get as limit 



2 



p x {v 2 )U(x, (y - v 2 ) + , k) dv 2 . 



From this procedure we pick up another o(e) term. 
The last case has probability 

r+oo /-+00 

/ / d y U(x,z,k- l)p x (vi)p x (y 2 )\{ Vl+V2 > y+z+£ }dvidv 2 dz. 

Je Je 

By assumption, —d y U(x,y,k — 1) exists and it is the density of the measure 
U(x,-,k). The dummy variables vi,v 2 , z stand for the excess heights of T_i, T\, Tq 
respectively, and in this case the central slope in A x+£ {w) has excess height V\ + 
v 2 — z — e, giving the restriction i>i + v 2 — z — e > y. Since d y U{x, -,k — 1) is 
continuous, dividing by e and taking e — > we get 

r+ca r+00 

-d y U(x,0,k-l) / l{ Vl +v 2 >y}Px{vi)p x (v 2 )dv 1 dv 2 . 

Jo Jo 

And again we pick up an o(e) term. The double integral equals (y/x+1) exp(—y/x). 



Before getting to the actual proof of our theorem, we give a non-rigorous short 
derivation to illustrate its main steps. The main problem is that we don't know if 
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U is differentiable in the x, y variables for every k G N. Assume for the moment 
that it is. For f,g G L 1 ([0, +00)), as usual, we define / * g G C([0, +00)) by 

(/ * g)(x) ■= Jo f( x - y)a{y) dy for ail x e [o, +00). 

The above lemma would give for U the PDE 

(14) (d x -d y )U(x,y,k) = -- U(x,y,k) + ^ (U(x, ■, k) * e~*)(y) 

x 

2 _y 

H — e " U(x, 0, k) — l{k>i}dyU{x, 0, k — l)(y/x + 1) exp(—y/x). 

Let /(a;, y, fc) = L/(x, j/x, fc); i.e., U(x, y, k) = f(x, y/x, k). 
Then / should satisfy 

(15) (x d x -(1 + y) d y + 2)f(x, y, k) = 2 (f(x, k) * e -)(y) + 2 e -*/(s, 0, k) 

-l {k > 1} (y + l)e- y d y f(x,0,k-l), 

while the conditions at x — 1 are 

(16) /(l,y,0) = (2y/3 + l)e- y for y > 0, 

(17) /(l,y,fc) = fory>0,fc>l. 

The first equation comes from 1|11[) . the second is clear. 

For z G D := {z G C : |z| < 1}, the generating function M(x,y,z) := 

oo 

f(x,y,k)z k is well defined. Assuming that M is differentiable with respect 

to x,y and its a;, y derivatives are obtained with term by term differentiation, we 
see that M satisfies the PDE problem 

(x d x -{l + y) d y + 2) M {x, y,z) = 2 (M (x, -,z)* e~)(y) 

(18) +2e~ y M{x,0,z)~(y + l)e~ v zd y M(x,0,z) in (1, oo) X (0, oo), 

(19) M(l,y,z) = (2j//3 + l)e~ v for y > 0. 
We try for a solution of the form 

M(x, y, z) — [a(x, z) + b(x, z) y] e~ y . 

Substituting this into l|18|l we see that e~ v factors out in both sides, and after 
cancellation, we arrive in an equality of two first degree polynomials in y with 
coefficients depending on x, z. Equating the coefficients in equal powers of y in the 
two sides of the equation, we arrive at the following system of ODEs for a, b. 

x d x a(x, z) + (z — 1) (b(x, z) — a(x, z)) = 0, 

x d x b(x, z) + (2 + z) b(x, z) — (z + 1) a(x, z) = 0. 

The initial condition 119|) for M, expressed in terms of a, 6, becomes a(l,z) = 
1, 6(1, z) = 2/3 for all z G D. We easily see that the only solution of the system 
satisfying these conditions is 

a(x, z) = ci(z) x Xl{z) + c 2 (z) x X2{z \ 

b(x, z) = Cl (z) + x^ + c 2 (z) (l + ^M) x^ z \ 

where 

. . . -3 + a/5 + 4z . . . -3 - a/5 + 4z 
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and 

Cl (z) = ((z-l)/3-A 2 (z))/(A 1 (z)-A 2 (z)), c 2 (z) = (-(z-l)/3+A 1 (z))/(A 1 (z)-A 2 (z)). 
Now E{z k ^) = M(x, 0, z) = a(x, z), which is what we want. 

Proof of the Theorem: The proof is done by taking the steps of the above 
"proof" in reverse order. This time all the steps can be justified. We will need 
three lemmata. Two of them are non-trivial, and their proof is given in section 

Lemma 3. The solution M of (|18|l . 119|) obtained above is analytic as a function 
of z in D. The coefficients of its development as a power series around zero are 
differ entiable with respect to x and y in (1, +oo) x (0, +oo), and its x, y-derivatives 
are continuous and can be found with term by term differentiation. 

Using this lemma, we write M as 

oo 

M(x,y,z) :=^2g{x,y,k)z k . 

fe=0 

Differentiating M term by term and equating the coefficients of equal powers of z 
in the two sides of (|18|) . we see that the sequence of functions [g(-, •, k)) k>Q satisfies 
the PDEs l(T5|l with conditions at x = 1 given by fTHjl , fT7|l . 

For k G N define g{-, •, k) : [1, +oo) x [0, +oo) — > K by g(x, y, k) — g(x, y/x, k) 
for (x,y,k) e [1, +oo) x [0, +oo). The sequence of functions •, fc)) fe>0 satisfies 
the PDEs with conditions at x = 1 given by fTH|) . fT7|l . 

The proof is finished by showing that the sequence (U(-, •, fc)) fceN satisfies a weak 
form of these PDEs, and then a uniqueness result will identify U as g. 

For fixed c > 1, define g c ,k(x) = U(x, c— x, k) for x £ [1, c]. We state as a lemma 
an equation that g Cy k satisfies. The proof is straightforward from Lemma 

Lemma 4. The function g c ^ is differentiable in (1, c) and satisfies 

(20) g' k (x) =-- g Cik (x) + \ (U(x, ■, k) * e"=) (c- x) 

X x z 

+ - e -^r U(x, 0, k) - lr fe>1> d v U(x, 0, k - 1) (cAr) e - {c - x)/x , 
x 

where for k > I we assume that U(x, y, k — 1) is differentiable in y with continuous 
derivative. 

And the promised uniqueness result is the following. 

Lemma 5 (Uniqueness). Let f : [l,+oo) x [0, +oo) — * [0, +oo) be a continuous 
function such that for every c > 1 the function g c : [l,c] — > M wii/i gc(a;) := 
f(x,c— x) is continuous on [l,c] a?i(i differentiable on (l,c) ; wii/i 

(21) ^(i) = -- ffe (i) + A (f(x, ■) * e--)(c ~x) + -e-^f(x, 0) 

X X z X 

and g c (l) = 0. Then f = 0. 

Now using induction we show that 

U(-,-,k) = g(-,-,k) VfceN. 

The function £/(■,-, 0) — g(-,-,0) satisfies the assumptions of Lemma because of 
the PDE problem that g(-,-,0) solves and Lemma 0J For k > 1, assuming the 



DIFFUSION IN RANDOM ENVIRONMENT 



15 



statement true for k— 1, the same argument works for {/(•, •, fc) — <?(•, •, fc), where now 
the assumption on {/(•,•, fc — 1) required by Lemma 01 is provided by the inductive 
hypothesis. 

Therefore, J2k=o xy, k) z k = M(x, y, z) for (x, y, z) e [1, +oo) x [0, +oo) x D 
and E(z k ^) = M (x, 0, z) = a(x, z), proving © for z e D. ■ 

Proof of Corollary |5J Our theorem gives that P(6 doesn't change sign in 
[l,x]) = a(x,0), where the function a is defined on pagc ll3l and by scaling 

P(6 doesn't change sign in [x, y]) = a(y/x, 0) for < x < y. 

Consequently, the density of the last point before y that we have sign change is 
—yx~ 2 d x a{y / x, 0), where we use d x here and below to denote derivative with respect 
to the first argument. Differentiating with respect to y, we get the density of the 
event that x, y are consecutive times of sign change as 

x~ 2 d x a(y/x, 0) + yx~ 3 'd xx a(y / 'x , 0), 

which, after using the expression for a(x,0), becomes 

The event that b changes sign at x translates to the central x-slope having excess 
height 0. This has density l/(3x) due to the scaling property of the x-slopes and 
relation Qll[l. which refers to the central 1-slopc. Thus, the density of the time Y 
where the next sign change after x happens, given that there was a sign change at 
x, is 



(22) h(y) 



otherwise. 



The quotient Y/x, given that there was a sign change at x, has the density of the 
fj's given in (JHJ. 

■ 

Proof of Corollary [1} Observe that the right-hand side of © is a function 
analytic in C \ (— oo, —5/4]. As for the left-hand side, we have the following lemma. 

Lemma 6. For any x > I, the power series 

F,n=o ¥ ( k (x) = n) z n defines an entire 

function. 

And from a basic property of analytic functions, it follows that the quantities 
M(x, 0, z), a(x, z) agree for all z € C \ (— oo, —5/4]. ■ 

Remark 4. Of course, Lemma El implies that the right hand side of Q can be 
extended to an entire function. However, the way © is written doesn't allow as to 
claim that it holds for all z e C because the function (z — > + 4z) doesn't have 
an entire extension. 

Proof of Corollary 151 Wc apply the Gartner-Ellis Theorem (Theorem 2.3.6 in 
0]). The moment generating function of k(e t )/t is given for any A € R by A t (A) := 
logE(exp{Afc(e t )/i})- So ^A^iA) = ^ logE(e AA; ( et )) = i" 1 logM(e', 0, e A ), and 
using Corollary 01 we see that 



1 a ,,„ , / an -3 + V5 + 4e A 



A(A) := lim -A 4 (tA) = Ai (e ) = 

t—^OO t 
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The Fenchel-Legendre transform A* of A, defined by A*(x) = sup AgK {Ax— A(A)} for 
all x G R, is found to be the function / defined in the statement of the proposition. 
Also Da '■= {A G R : A(A) < 00} = R, and A is strictly convex and differentiable in 
Da- The result follows from the Gartner-Ellis theorem. ■ 

4. Proofs of the lemmata 
Lemma 7. For any x > 1, y > 0, k G Z, and e > 0, we have 
\U(x + e, y, k) — U(x, y + e, k)\ < 3e/x, 
\U(x, y, k) — U(x, y + s, k)\ < e/x. 
In particular, U is continuous. 

Proof. Call A and B the two events whose probabilities are U(x+e, y, k) and U (x, y+ 
e,k) respectively. Then AAB C [ In A x (w) at least one of the three slopes neigh- 
boring zero has excess height < e ]. Denote by Tq,T\ the central slope and the 
slope to the right to it in A x (w). Then 

F(AAB) < 2P( Ti has excess height < e ) + P( T has excess height < e ) 

ne/x re/x 

= 2 1 e- z dz+ (2z/3+ l/3)e" 2 dz < 3e/x. 
Jo Jo 

The other inequality follows because the density of the measure lA{x, S, k) is bounded 
by 1/x (sec ■ 

Proof of Lemma|3j Define K : CxCxD -> C within (zi, z 2 , z) := M(e Zl ,z 2 ,z) 
for (zi, Z2,z) £ CxCxfl. Clearly K is a holomorphic function (choose an analytic 
branch of the square root function defined on C \ (—00, 0); the number 5 + 4z is 
there for z G D) and has a power series development centered at zero that converges 
inCxCxD (see, e.g., |2J, Proposition 2.3.16). The claims of the lemma follow by 
the relation M{x, y, z) = K(logx, y, z) and standard properties of power series. I 

Proof of Lemma[5j For c > I and x G [1, c], define N(c, x) := sup{|/(z, c—z)\ : 
z G [1,2;]}. N is well defined because / is bounded on compact sets. 
From (|2f \ , for x G ( I , c) one has 

\g' c (x)\ < 2 N(c,x) +4 sup N{d,x), 

l<d<c 

and since g c (l) = 0, we get after integrating 



\f(x,c-x)\ <6 



/ sup N(d,t)dt, 

Jl Kd<c 



which implies 



N(c,x) < 6 / sup N(d,t)dt 

'l l<d<c 



and 

sup N(d,x)<6 / sup N(d,t)dt. 

l<d<c Jl l<d<c 

The function A(x) = sup 1<d<c N(d, x) (x G [l,c]) is continuous (because / is) and 
has A(l) = 0. An application of Gronwall's lemma to A gives N(d,x) = for 
x G [I, c], d G [I, c]; i.e., f(x, y) = for all (x, y) G [I, +00) x [0, +00). ■ 
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Proof of Lemma |5J Let Xk be as in Corollary Then 

P(fc(x) = n + 1) < P(fc(x) > n + 1) = P(X n+J < ar) 

< P(rir 2 ■■■r n <x)= P(logn + log r 2 H h logr„ < logo;). 

For i > 1 set Yj = log r,, 5; = Yi + I2 + ■ ■ • + Yi , and let be the dis- 
tribution measure of By Cramer's theorem (0] Theorem 2.2.3), the se- 
quence (mji>i satisfies a large deviation principle with rate function I(x) = 

sup AeR {Aai -TogE(e Alogri )} for x G R - 
Clearly, for any £ > we have 

IW ilogP(^ < ^) < -7(e). 

n— foo n n n 

And lim e _>.o -f(e) = +°o because I is lower semicontinuous and 1(0) = +00. To see 
the last point, observe that 

7(0) = sup{-logE(e Alogri )} > BE {-logE(e Alogri )} = +00 

because limA_>_oo e Alogri = with probability one and the bounded convergence 
theorem applies. Thus Unin—xx, — logP(fc(a;) = n + 1) = —00, proving that the 
radius of convergence for the power series is infinite. ■ 

Lemma 8. For P and Wi as defined in the introduction, it holds that P(VVi) = 1. 

Proof. First we prove that for fixed x > 0, the set 

C x := {w 6 C(R) : R x {w) has the properties appearing in the definition of Wi} 

has ¥(C X ) = 1. Observe that for z a point of x-minimum and a z := sup{a < 
z : w(a) > w(z) + x}, Z := inf{/3 > z : w(f3) > w(z) + x} it holds that a z < 
z < f3 z because w is continuous at z. And there is no other point of x- minimum 
in (a z ,/3 z ) since if z is such a point, say in (a 2 ,z), then assuming /3g < z we get 
a contradiction with the definition of a z while assuming f3 z > z we get that w 
takes the same value in two local minima, which has probability zero. Now assume 
that there is a strictly monotone, say increasing, sequence (z n ) n >i of x-minima 
converging to z^ G R. By the above observations we get lim yj yy Zoo (w(y) — w(y)) > 
x contradicting the continuity of w at z^. Similarly if (z„)„>i is decreasing. So 
in a set of w's in C(K) with probability 1, it holds that the set of x-minima of 
w has no accumulation point. The same holds for the set of x- maxima and as a 
result also for R x (w). Since lirri | f | >OQ w t = —00, lbxi\t\^oo w t = +00, it follows 
that P(R x (w) is unbounded above and below) = 1. ft is clear that between two 
consecutive x-maxima (rcsp. minima) there is exactly one x-minimum (rcsp. x- 
maximum). Consequently, P(C X ) = f . 

Finally, note that for all n £ N \ {0} we have R n (w) C R x (w) C Ri/ n (w) for 
x £ [l/n,n], from which it follows that Wi = n xe (o )+00 )Cx = n, lS N\{o}( C '» nC 'i/n)- 
Thus, P(Wi) = f . ■ 
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